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ABSTRACT 

Bayesian statistical methods of Gregory-Loredo and the Bretthorst generalization of 
the Lomb-Scargle periodogram have been applied for studying activity cycles of the 
early K-type subgiant star CF Oct. We have used a ~ 45 year long dataset derived from 
archival photographic observations, published photoelectric photometry, Hipparcos 
data series and All Sky Automated Survey archive. We have confirmed the already 
known rotational period for the star of 20.16 d and have shown evidences that it 
has exhibited changes from 19.90 d to 20.45 d. This is an indication for stellar surface 
differential rotation. The Bayesian magnitude and time-residual analysis reveals clearly 
at least one long-term cycle. The cycle lenght's posterior distributions appear to be 
multimodal with a pronounced peak at a period of 7.1 y with FWHM of 54 d for 
time-residuals and at a period of 9.8 y with FWHM of 184 d for magitude data. 
These results are consistent with the previously postulated cycle of 9 ± 3 years. 

Key words: methods: statistical - stars: individual: CF Oct - stars: activity - 
astronomical databases: miscellaneous. 



1 INTRODUCTION 

Long-term stellar activity study is a challenging research 
area so far as it needs the longest available data sets, 
and appropriate statis tical methods to treat the problem- 
atic in its complexity. I Phi lips fc HartmannI (|l978l ') first re- 
ported the long-term photometric variability of spotted bi- 
nary stars, based on the archival photographic plate photom- 
etry combined with photoelectric photometric data. Their 
work produced more than 50 years of activity history for 
BY Dra and CC Eri. A sim i lar re sult for V833 Tau was 
pu blished b y Hartmann et al.l ([l98l|), and is later confirmed 
bv lBondaJ (| 19951 ). Star-spot activity is now known to oc- 
cur on a range of solar-type and early K giants stars. Some 
releva n t works i n the field are those of [ Messi na fc Guina n 
ll2002h.lB ondar' ('l995 M01ah. Kollath fc Strassmeieij (|2000i ). 
lOlah fc Strassmeier, (,20o3 ). 

In the present analysis we used combined photographic, 
photoelectric photometry and CCD observations to study 
cyclic activity of the spotted subgiant star CF Oct. The 
collected data, presented on Fig. 1, are quasi-randomly dis- 
tributed with data gaps over time interval of about 45 years. 
We applied the Gregory-Loredo Bayesi an method for t ime 
series with independent Gaussian noise, ICregorvl (|l999l ). m 
order to obtain Bayesian estimate for the known ~ 20 d 
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rotational modulation and to search for long-term periodic 
variability in time-scales up to about 15 years. The method 
and the IDL procedure, developed in the AIU Jena, that we 
used have been previously tested on the synthetic data-sets, 
modeled with known period and amplitude, observational 
errors and simulated data gaps. They have shown reliable 
results and accurate period estimation. 



2 CF OCT - ACTIVE SPOTTED GIANT STAR 

CF Oct (HD 196818, HIP 102803, a = 20''49^2371263jJ = 
-80°0 8'01".01, J2000, KOIIIp given by iHouk fc Co^devl 
(| 19751 )) is a bright, southe rn active gi a nt sta r. Its variabil- 
ity was first detected by IStrohmeieJ (| 19671 ) on the pho- 
tographic plates from Bamberg Observatory Sout hern Sky 
Survey (BOSSS). In the GCVS (|Samus et al.ll2010l ) the star 
is described as a RS CVn type variable with maximal bright- 
ness V = 8.27 mag and relatively large photometric vari- 
ations ~ 0.3 mag. Early archival observation s from the 
BOSS S were recently digitized and analyzed bv llnnis et al.l 
(|2004 ). 

Photoelectric photometry for t h e star was performed 
lm§i]. iLlovd Evans fc Koen 



Innis. Coates fc Thompson 



and pub lished by [ Innis et al 
(198^, iPollard et al.l (| 19891 ) . 
(19971 ). Rap id changes in the ligh t curv e in 2006 was an- 
nounced by llnnis. Coates fc Kavd (|2007l ) . The photometry 
studies reveal rotational modulation of ~ 20 d, due to spot- 
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ted activity . Spec troscopic studies bv lHearnshawl (|l979l ) and 
llnnis et al.l l| 19971) of CF Oct s how s trong Ca ii emission and 
filled-in Ha line, llnnis et all (|l997l ) reported for no signif- 
icant radial velocity variations and hence concluded that 
the star is probably single. Based on the radial velocity 
data and the period of the rotational cycle they consider 
it is a subgiant rather than a giant as per Houk classifica- 
tion. CF Oct is also re ported to b e a st rong, flaring, mi- 
crowave radio source by ISlee et all lll987h. and appears at 
the ROSAT Bright survey catalogue~ (|Fischer et al] Il998l ') 
with 1.12 counts per second in the energy range of 0.1- 
2.4 keV. 



3 OBSERVATIONS 

For the present study we have used data from: the Bam- 
berg Observatory Southern Sky Survey (BOSSS), published 
photoelectric photometry observations, the Hipparcos satel- 
lite time-series, and the All Sky Automated Survey (ASAS) 
data archive. 

BOSSS l|Tsvetkov et al.ll2005h was taken in the period 
1963 to 1976 and contains more than 22 000 plates that 
cover the whole southern sky with limiting magnitude in 
range 11 to 14 mag. This survey is unique because of the 
fact that in this time interval the first Palomar and Har- 
vard Observatory Southern Sky Surveys were completed, 
but other southern observatories were not actively surveying 
the sky. C F Oct was observe d on 352 plates in the period 
1964-1976. llnnis et al.l (|2004l ) performed aperture photom- 
etry and p hase dispers i on mi nimization (PDM) light curve 
analysis. In llnnis et all (|2004l ) we discussed the methods for 
photographic plate digitization, the plate photometry and 
the transformation from B toV magnitude of the star. Plate 
data (BAM) are collected in more then six seasons relatively 
well covered by observations, separated by significant inter- 
vals with a lack of observations. 

Photoelectric photometr y data (PHOT) came from 
follow i ng papers: j lnnis ct al.' (1983'), 'Llovd E vans fc KoenI 
l|l987h . |Pollard et al., (1989) , Innis et al , (1997). The dataset 
consists of 137 data points, that might be separated in 
5 seasons of intense observations. In this PHOT dataset 
the mean brightness of the star is at its historical min- 
imum. It is no t clear whether this is an observational 
selection effect. iLlovd Evans fc KoenI l| 19871 ) applied sys- 
tematic magnitude corrections in order to have agreement 
for the photometry from different seasons, although they 
noted this was purely an ad hoc approach. Changes of 
mean light level have been observe d for other RS CVn 
stars (iFrohlich. KroU fc Strassmeierl l|2006l ) for HK Lac, 
iPoretti et all l| 19861) for AY Cet). 

The Hipparcos and Tycho Catalogues IESAI (|l997h . 
available via the Centre de Donnes astronomiques de Stras- 
bourg (CDS), provide high quality scientific data for the pe- 
riod November 1989-March 1993. For our research we have 
used the Hipparcos time series data, as they have relatively 
smaller errors and are close to the photometry from Ty- 
cho catalogues. The Hipparcos magnitudes of CF O ct were 
transf ormed to V by the use of equations given by iBessell 
l|2000l ). 

The most r ecent photometry for our work is extracted 
from the ASAS. IPoimanskH (|2002l ). archive at the web-page 






1.0x10^ 
HJD - 2436000 



Figure 1. V magnitude liglit curve of CF Oct for the period 
1964—2009 with overploted errors in gray colour. 



Table 1. Statistical description of the data 



Dataset 


Np 


Ts 


HJD 




Vmax 


<V > 


BAM 


352 


4484 


2438560 


7.67 


8.30 


7.98 


PHOT 


137 


3212 


2444071 


7.93 


8.41 


8.16 


HIP 


130 


1176 


2447873 


7.74 


7.98 


7.90 


ASAS 


705 


3058 


2451961 


7.73 


8.07 


7.91 



http: / /www. astrouw.edu.pl/asas/ The ASAS project aims 
at photometric monitoring of the whole available sky, from 
two observing stations - Las Campanas Observatory, Chile 
since 1997 and Haleakala, Maui since 2006. Both stations are 
equipped with two wide-field 200-mm/f/2.8 instruments, 
observing simultaneously in V and I band. The project's 
ultimate goal is detection and investigation of any kind of 
photometric variability. As CF Oct is a relatively bright star, 
so far the early ASAS photometry data are unusable due 
to saturated images. For this reason, after consulting with 
ASAS supporting team, for the analysis we included only 
the data free of saturation infiuence. 

The dataset we have collected contains data for HJD of 
the observation, V mag of the star, and corresponding errors. 
As far as it consists of four datasets taken from different 
sources it suffers from significant intervals with lack of data, 
and also by non-uniform data distribution. The resulting V 
magnitude light curve with overploted errors is presented 
on Fig. 1, where the crosses present BAM data, asterisk - 
PHOT data, diamonds - HIP data and triangles - ASAS 
data. The data-set statistically is presented in Table 1, with 
following information: Dataset ; A^p - the Number of data 
points in the set; Ts - the time span of the set in days; HJD is 
in the beginning of the set; Vmin, Vmax and <V>- minimal, 
maximal and mean values of V magnitudes respectively. 



BAYESIAN STATISTICAL METHODS IN 
USE 



We used Gregory-Loredo methodfCL). ICregorvl (|l999l ). for 
time-series analysis in the presence of Gaussian errors. It 
employs Bayes's theorem for detection of a periodic sig- 
nal with unknown shape and period and for estimation of 
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the characteristics of the detected signal. First step of the 
Bayesian statistical methods is to determine the hypothesis 
space. In general, for observed stars we have three hypothe- 
ses - constant light, variable light, and a periodic signal. In 
the GL method periodic models are represented by a signal 
folded into a stepwise function, similar to a histogram, with 
m phase bins per period plus a noise contribution. With 
such a model we are able to approximate a light curve of 
any shape. 

The stepwise model is then described by the following 
parameters: period P, phase of the minimal star brightness 
(j), number of bins m, light curve shape ri in the each bin 
and noise scale parameter b. The GL method uses Jeffreys 
prior distributions for b and P, to ensure fully compatibility 
of the posterior distribution functions of the period and the 
frequency (1/-P) , and uniform prior distribution for the rest 
of the parameters. The prior range for the period and the 
light curve shape parameter are user selected according the 
available data and the prior information. 

The constant model may be treated just as a special 
case of a periodic model with only one bin m — 1 and it has 
only two parameters: the value of the constant signal and 
the noise scale parameter. In the stepwise representation the 
nonperiodic modulation model may be treated as another 
special case of periodic model, when the period is equal to 
the data time span, with wider range of the number of bins 
and it has m light-curve shape parameters, bin phase and 
the noise scale parameter. 

The stepwise model is flexible enough to approximate 
light curve of practically any sha pe, and althoug h it is well 
suited to describe "spiky" signals, lGregorvll|l999l ) points out 
that it can successfully detect smooth signals and accurately 
estimates their shape and parameters. The posterior mean 
values, the most probable values and the credible intervals of 
the parameters are then estimated by marginalization of the 
global likelihoods over the priori specified range of each pa- 
rameter. In our particular case the posterior mean is a robust 
estimator for the rotational period the of CF Oct. Bayesian 
posterior probability contains a term that penalizes com- 
plex models, hence we calculate the posterior probability by 
marginalizing over a range of models, corresponding to a 
prior range of m from 2 to 12. The noise scale parameter 
b in GL method is the relation of the variances of the ob- 
server's errors and model residuals, and it is an indication 
whether the model successfully describes the observational 
data. If the mode of the 6 ~ 1, then the proposed model 
is accounting for everything that is not noise in the data. 
If the mode of the & < 1, then the model does not account 
significant features in the data, or the initial observers error 
estimates are low. Values of the & > 1 is an indication for 
over-fitting as the model residuals are smaller than the ob- 
servational errors. We have tested the method by using mod- 
eled datasets, with randomly distributed data-points and 
with data gaps. For independent verification of the detected 
signals we have also applied Lomb-Scargle method, together 
with the Ba yesian gene r alizat ion of the Lomb-Scargle peri- 
odogram bv lBretthors^ (|200ll ). 

In order to reduce the likelihood of introducing un- 
known systematic errors into the periodicity study of 
CF Oct, due to the different observational methods and 
data reduction procedures employed, instead of analyzing 
V magnitude variations, we have analyzed the variations of 



< V > —V, where < 1/ > is the mean value of V magni- 
tude for each data set, i.e. < V > = (< Vbam >,< Vphot >,< 
Vhip >,< Vasas >)• On a short time-scale the star is already 
known as a variable, so while we analysed the rotational 
modulation we concentrate on testing the hypothesis for de- 
tecting periodic variations, and on the period estimation. 
With the collected data we are able to study variability of 
CF Oct in ranges from several days (this limit is set up by 
the average sampling frequency of our observations) to 15 
years (1/3 of the covered observational time span). This is 
a rather large time span, and was examined in several parts 
by using a suitable number of frequencies in each section. 
The long-term variability study of the star was performed 
by estimating the odds ratios (the ratio of the probabili- 
ties of the model with respect to the constant model) of the 
nonperiodic and periodic modulation models, as well as by 
searching for periodic modulation. 



5 DETECTING THE ROTATIONAL 
MODULATION 

As CF Oct has a well established rotational period near ~ 
20 d, known from PDM and least-squares periodical analy- 
sis of photographic and photoelectric observations, we have 
applied the Gregory-Loredo method for its detection. This is 
relevant as already studie d BAM and PHOT data sets show 
slightly different periods (|lnnis et al.ll200j ') that might be 
due to the period estimation errors, phase shifting or more 
complex periodical modulation. All the data sets are rich 
enough and are suitable to be studied separately as well as 
all together. For the separate datasets and for all the com- 
bined data we have calculated the global posterior probabil- 
ity density functions for a class of models described by the 
following parameters : 6 noise scale in the priori range from 
0.05 to 2.0, number of bins (m) in range from 2 to 12; the 
period (P) in a prior range from 19 to 21 d, cove red with 
relevant number of frequencies. We have followed ICregorvl 
(1999") in the case of fe 7^ 1 and estimated the projected prob- 
ability density function of the noise scale parameter. In case 
of the PHOT and ASAS data we have used a model with 
additional error term, in order to avoid numerical difficulties 
because of the very low observational errors. 

The parameters of interest are estimated by margi nal- 
ization of the global likelihood function (equation 7. iGregorvl 
(|l999l) ) over the nuisance parameters in the priori range. 
Marginalization over the noise scale parameter is the most 
conservative way for estimation of the model parameters 
that treats anything in the data, not described by the model 
as a noise. By marginalization of the joint posterior likeli- 
hood function over the frequencies we have computed the 
probability of the number of bins, and thus the most prob- 
able m (mmodo in Table 2). The number of bins parameter 
relates to the complexity of the light curve to the shape of 
the light-curve and probably is connected with the struc- 
ture of the stellar spots (or the spot groups). The values of 
the mmodc parameter, we obtained, tend to be small 2 or 3 
with a rather narrow distribution, so that model with one 
more light curve bin is already penalized by the Occam's ra- 
zor factor for involving unreasonable number of parameters. 
The marginalization of the joint probability function over 
m represents the information provided by the data for the 
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Figure 2. Marginal PDFs (see the text) of the period, plots show only the interesting part {left panel) and mean V subtracted light 
curve , with common HJD = 2438560.4 for the zero phase {right panel) for: from top to bottom - BAM, PHOT, HIP, ASAS and ALL 
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Table 2. Derived parameters: number of bins, period and 68 per cent 
credible interval 



Dataset 


^modc 


^modc 


-^modc 


p 

^ mean 


ere. int., (from — to) 


^brct 




Value 


(days) 


(days) 


(days) 


(days) 


(days) 


BAM 


3 


0.95 


20.04 


20.04 


20.035-20.045 


20.04 


PHOT 


2 


0.80 


20.17 


20.16 


20.145-20.175 


20.16 


HIP 


2 


1.00 


20.46 


20.45 


20.425-20.475 


20.51 


ASAS 


3 


0.90 


19.94 


19.94 


19.930-19.950 


19.93 


ALL 


3 


0.28 


20.16 


20.14 


20.120-20.950 


20.11 



period. The period that best describes the observed hght 
curve may be estimated by determining the maximum of 
the marginal posterior probabihty dencity function, or by 
calculating the posterior weighted mean period. We have 
calculated the 68 per cent credible interval s for a given pa- 
rameter following ICregorv fc Loredd (jl992j ) as the interval, 
that contains 68 per cent of the posterior Probability Density 
Function (PDF) and where the PDF is everywhere greater 
than the one outside the credible interval. The Bretthorst 
generalization of the Lomb-Scargle periodogram gives very 
close results and also confirms the period variations. Table 
2 represents our results: the most probable mmodc, the most 
probable noise scale parameter &mode, the maximal prob- 
able period (Pmode), the weighted mean period Pmoan, the 
credible intervals for the period detection and the period de- 
rived by Bretthorst generalisation method Pbrct ■ The poste- 
rior probabilities are normalized in order to obtain the PDFs 
and the later have an integral over the priori ranges equal 
to 1. For the separate datasets as well as for all the data 
together the expanded plots showing the interesting part of 
the marginal PDFs are given on the left panels of Fig. 2. The 
mean subtracted V light curves, plotted with the most prob- 
able period for the datasets and for all the data, and with 
a common HJD for the zero phase for all the photometries, 
set at the beginning of observations at HJD = 2438560.4 
are presented in the right panels of Fig. 2. For the separate 
datasets ^modc SuS well 3.S tll6 &mcan are close to 1, and shows 
that the periodic model is doing a good representation of the 
obser vational data. Ou r analysis confirms previous sugges- 
tions (jlnnis et al.ll2004 ) for period change s as well as for light 
curve phase shifting TPollard et al.lll989l ). The value of the 
femodc ~ 0.28 for ALL the data indicates that the simple pe- 
riodic model does not explain satisfying all the observational 
data and eventually refers for the presence of additional sig- 
nal in the data. However, in the case of combined data we 
do see possible other source of "noise" (on the base of the 
estimated parameter b) and in principle larger credible inter- 
val that depends not only from spanned time, but also from 
the signal to noise ratio and from the amplitude of periodic 
variations. The mean and mode period obtained for photo- 
graphic plate and photoelectric photometry data are close 
to the previous published ones. Scatter in the light curves 
from different datasets shows that the amplitude of bright- 
ness variations, light curve shape and phase of the minimal 
brightness change with the epoch of observation. This addi- 
tionally complicates the estimation of the light curve shape 
which is beyond the goals of the presented research. 



6 LONG-TERM VARIABILITY OF CF OCT 

As it is mentioned above, the GL method gives us oppor- 
tunity for searching for long-term cycles with duration up 
to 6000 days. For the proposed study we selected time-scale 
from 600 to 6000 days that is rather different from the time- 
scale of the rotational modulation and allows to separately 
study the long-term variability of CF Oct. The hypotheses 
for constant (He), for nonperiodic (//np) and for periodic 
signal (Hp) are evaluated, and the odds (hi/he) ratios are 
given in Table 3. The noise scale parameter for all the mod- 
els takes into account the rotational modulation as an ad- 
ditional noise. For the nonperiodic model we use number of 
bins parameter m — 40. The result shows that the periodic 
model can best represent the observational data and it is 
the most probable model. The variable model appears the 
second probable one and it is very reasonable as periodical 
models are a special case of variable ones. The noise scale 
parameters for the examined models give the additional in- 
formation that the constant and variable models on a long- 
term scale do not explain significant structures in the data. 

In order to search for a long-term periodicity we have 
computed the Bayesian PDF and estimated probabilities for 
period detection over the mean-subtracted magnitude (< 
V > —V , where < V > is the mean value of V for each data 

set, i.e. < V > = (< Vbam >,< Vphot >,< Vhip >,< Vasas >) 

data on long-term time scales with a prior period range from 
600 to 6000 d. Bayesian estimation of the probability for 
long-term periodic modulation over mean-subtracted mag- 
nitude data results in a marginal PDF (Fig. 3) with three 
individual peaks with local maxima at 3582 d (~ 9.8 yr), 
2432.5 d (- 6.7 yr) and 1173 d (- 3.2 yr) and FWHM of 
184 d, 70 d and 21 d respectively. 

The phase shift of t he mi nimal or maxima l light sug- 
gested by iPollard et all l|l989l ) and llnnis et al.l (|2004l ') has 
led us to the idea to study time-residuals. The time-residuals 
represent the difference between the observed and model 
predicted times of signal minima (or maxima) and are an 
analog of the (O — C) terms in light curves of the binary 
stars. They are representing the observed phase shift in the 
light curve of CF Oct and are probably connected with the 
spot migration over the stellar equator. We have computed 
the time-residuals by using stepwise (with m = 10 bins) light 
curves of the PHOT dataset as a reference for expected times 
of minimal brightness of the star. The differences between 
the expected or predicted by the PHOT light curve and 
actually observed times of minimum brightness for all the 
observed data span were calculated. The marginalized (over 
m) probabilities for detection of cyclic activity with period 
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Table 3. The hypotheses space, odds ratios and noise scale parameters 



Hypothesis 


odds ratio hi/ho 


^modc 


He : Constant signal with randon Gaussian noise 

H]\fp : Nonperiodic signal with unknown shape 

Hp : Periodic signal with a priori unknown shape and period 


1 

1.8 X 10^ 
0.4 X 10**" 


0.21 
0.54 
1.32 



Figure 3. Marginal PDF of the period for mean-subtracted V 
magnitude data 



Figure 4. Marginal PDF of the period for time-residual data 



in the interval 600 to 6000 d are presented on Fig. 4. The 
marginal total probability of 0.46 for cycle with a period of 
2603.0 d (~ 7.1 yr) with FWHM of 54 d is estimated. This 
cycle length i s in a good agreement with the suggestions of 
IPoUard et al] (|l989l '). A period of 2648 d was detected with 
the Bretthorst generalization method and falls very close 
to the period determined by GL method. Fitting the time- 
residual data with a stepwise curve with period of 2603.0 d 
is presented on Fig. 6. 

The other two peaks of the marginal PDF with peri- 
ods near 4100 and 5500 d and FWHM of 432 and 769 d 
have relatively low probabilities both of 0.27 and they give 
much more scatter of the stepwise fitting curve compared to 
the 2603.0 d period. We suppose that the minor peaks are 
probably due to the data gaps and the quasi-random data 
distribution. Thus we consider that the period of 2603.0 d 
is the most probable for the time-residual data. 



7 DISCUSSION 

The results from Bayesian detection of the rotational mod- 
ulation for CF Oct are in good agreement with previously 
published values. Using the Bayesian statistical method 
meant we were able to estimate the probabilities for rota- 
tional modulation over the priori selected range, based on 
the previous studies, and to determine the mode and mean 
periods derived from separate datasets as well as for all the 
observational data. The rather narrow credible intervals for 
period detection over the four separate datasets, state the 
range where the simple periodic model best describes the 
observations and thus we propose that the star probably 
exhibit variations in the rotational period. We have tested 
the hypotheses for constant, simple periodic and variable 
model for all the observed data and find out that the pe- 
riodic model is still the most probable one. Obviously this 
model does not explain the changes in the rotational pe- 
riod, so more complicated analysis with expanded hypothe- 
ses space to account for both rotational and long-term ac- 
tivity of CF Oct need to be performed. We were able also to 
estimate the mode period of 20.16 d for all the collected data, 
but with relatively large credible interval. This is very close 
to the best pe r iod es timated over electrophotometric data by 
iPoUard et all(|l989h . The estimated rotational period varies 
from 19.94 d for ASAS data to 20.46 d for HIP data. The 
most evident hypothesis for changes in the rotational period 
is differential rotation of the stellar surface, and latitude mi- 
gration of surface spots. It is noteworthy to point out that 
ra pid changes in the fight curve of CF Oct were observed 
bv llnnis et all (|2007l ) in 2006. ASAS data set does not have 
enough observations at the time of these changes that would 
give additional information. However for the ASAS observa- 
tions after the reported changes near HJD = 2454040.0 the 
am plitude of the variatio n is relatively small, about 0.2 mag. 

iPoUard et all (|l989i) suggested for activity cycle of 9 ± 3 
years based on the changes of the amplitude of the spot- 
wave. The period for cyclic activity of 2603 d estimated with 
GL method by the time-residual data is in an agreement 
with this prediction. This new found cycle is another indica- 
tion for changes in the rotational period of the star. Inspec- 
tion of the stepwise fight curve with period 2603 d and (Fig. 
6) overploted on the time-residual data shows, that charac- 
teristic features - minima and maxima of the hght curve are 
well supported by the observations. The derived period for 
cyclic activity of 2432.5 d over the mean-subtracted magni- 
tude data is close to the period derived from time-residual 
data, but the FWHMs of these cycles does not overlap. 
Thus the hypothesis of the existence of the joint period for 
the magnitude and time-residual cycles in the presence of 
rotational modulation is very interesti ng and w i ll be tested 
later by the use of the method used bv I Gregory! l| 19991 ). The 
estimated cycle of 3582 d (Fig. 5) by the use of magnitude 
data according the represented parameters is the most prob- 
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Figure 5. Stepwise, phase-sliifted liglit curve for the cycle with 
period of 3582 d plotted over the mean-subtracted V magnitude 
data in gray colour. Symbols are the same as on Fig. 1. 




Figure 6. Stepwise time-residual fitting curve for the cycle with 
period of 2603 d overploted over the time-residual (Trez) data. 



CF Oct. Two cycles of 2603 d and 3582 d with FWHMs 
of 53 and 184 d were detected by the use of two different 
statistical techniques. First one cycle of 2603 d represent 
the changes of stellar magnitude, while the second one is 
related to the changes of the phase of minimal brightness of 
the star, so it means that there is no superposition of the 
two cycles. 

We have given arguments that support previous sugges- 
tions for smooth period changes, due to differential rotation 
and possible spot migration, and derived the most probable 
periods for different intervals of the hght history of CF Oct. 
The estimated rotational period are: 20.04 d for the interval 
1964-76 from Bamberg observatory plate data, 20.16 d for 
the interval 1979-88 from electrophotometry data, 20.45 d 
for the interval 1989-93 from Hipparcos data and 19.94 d 
for the interval 2000-09 from ASAS archive. The Bayesian 
analysis indicates that the hypotheses for presence complex 
cyclic activity of CF Oct is also relevant and it is a subject 
of future work to evaluate the proposition for observation of 
several long-period cycles over the joint analysis of magni- 
tude and time-residual data. 

This work demonstrates the value of the astronomical 
data archives and draws attention to preserving archival ob- 
servations for future exploitation. With the extension of the 
observational time series by the photographic plate data we 
are able to study decadal variations of bright active stars. 
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able one, and also f its into the suggested cyclic activity by 
iPoUard et all j 19891 ). 

The logio(Prot/Pcyc) = "2.1 index for CF Oct, for 
Pcyc = 2603 d, is in a good agreement with the values for 
active stars summarised bv lOlah et al.l l|2000l ). nevertheless 
on fig. 13 from the mentioned paper CF Oct would lie very 
close to the similar RS CVn star HK Lac an d would be also 
close t o the line of the main sequence stars of* Baliunas et all 
l|l996h . On the basis of other detections of long-term vari- 
ati ons with two to three cycles for several stars reported 
( 200C ) , and for the activity cycles analysis of 
( 20061 ) for HK Lac, we also propose that the 



bvlOlah et al 



iFrohlich et al 



hypotheses for complex long-term variability with superpo- 
sition of two or more cycles of CF Oct are relevant and needs 
future evaluating by the means of Bayesian statistics. 



8 CONCLUSIONS 

The main result of our work is that it clearly shows the 
evidence of a long-term activity cycle of the active giant star 



REFERENCES 

Bahunas S. L., Nesme-Ribes E., Sokoloff D., Soon W. H., 

1996, ApJ, 460, 848 
Bessel M., 2000, PASP, 112, 961 
Bondar N. L, 1995, A&AS, 111, 259B 

Bretthorst G. Larry, 2001, "Generalizing the Lomb- 
Scargle periodogram" , in Bayesian Inference and Maxi- 
mum entropy methods in Science and Engineering, Ali- 
Mohammad-Djafari Ed., pp 241-256, AIP conference pro- 
ceedings 

ESA, 1997, The Hipparcos and Tycho catalogues, ESA SP- 
1200 

Fischer J.-U., Hasinger G., Schwope A. D., Brunner H., 
Boiler T., Triimper J., Voges W., Neizvestnyj S., 1998, 
AN, 319, 347F 

Frohlich H.-E., KroU P., Strassmeier K. G., 2006, A&A, 
454, 295 

Gregory P. C, Loredo T. J., 1992, ApJ, 398, 146G 
Gregory P. C, 1999, ApJ, 520, 361 

Hartmann L., Bopp B. W., Dussault M., Noah P. V., 
Klimke A., 1981, ApJ, 380, L85 



8 A. P. Borisova, V. V. Hamharyan and J. L. Innis 



Hearnshaw J. B., 1979, Proc. lAU CoUoq. No. 46, eds, 
Bateson, F. M., Smak J., Urch I. H., p.371 

Houk N., Cowley A. P., 1975, Michigan Catalogue of Two- 
Dimensional Spectral Types for the HD stars. Vol. 1, (Uni- 
versity of Michigan, Ann Arbor) 

Innis J. L., Coates D. W., Dieters S. W. B., Moon T. T., 
Thompson K., 1983, IBVS, No. 2386 

Innis J. L., Coates D. W., Thompson K., 1997, MNRAS, 
289, 515 

Innis J. L., Borisova A. P., Coates D. W., Tsvetkov M. K., 

2004, MNRAS, 355, 591 
Innis J. L., Coates D. W., Kaye T. J., 2007, IBVS, No. 5762 
Lloyd Evans T., Koen M. C. C, 1987, SAAO Circ, 11, 21 
Messina S., Guinan E. F., 2002, A&A, 393, 225 
Olah K., KoUath Z., Strassmeier K. C, 2000, A&A, 356, 

643 

Olah K., Strassmeier K. C, 2002, AN, 323, 361 
Philips M. J., Hartmann L., 1978, ApJ, 224, 182 
Pojmanski C, 2002, Acta Astronomica, 52, 397 
Pollard M. J., Hearnshaw J. B., Gilmore A. C, Kilmartin 

P. M., 1989, JAp& A, 10, 139 
Poretti E. et al., 1986, A&A, 157, 1 

Samus N. N., Durlevich, O. V et al., 
2010, GCVS database, Version 20 10 Jan, 
http://cdsarc.u-strasbg.fr/viz-bin/cat7B/gcvsl 

Slee O. BV,"Nerson G. I," Stewart'^irT.7Wright" A. E., Innis 
J. L., Ryan S. G., Vaughan A. E., 1987, MNRAS, 229, 659 

Strohmeier W., 1967, IBVS, No. 178 

Tsvetkov M. et al, 2005, PASRB, 5, 303T 

This paper has been typeset from a T^rjX/ I^Ti^X file prepared 
by the author. 



